An anchoring transition at surfaces with grafted liquid-crystalline chain molecules 
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The anchoring of nematic liquid crystals on surfaces with grafted liquid crystalline chain molecules 
is studied by computer simulations and within a mean field approach. The computer simulations 
show that a swollen layer of collectively tilted chains may induce untilted homeotropic (perpendic- 
ular) alignment in the nematic fluid. The results can be understood within a simple theoretical 
model. The anchoring on a layer of mutually attractive chains is determined by the structure of 
the interface between the layer of chain molecules and the solvent. It is controlled by an interplay 
between the attractive chain interactions, the translational entropy of the solvent and its elastic- 
ity. A second order anchoring transition driven by the grafting density from tilted to homeotropic 
alignment is predicted. 
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Nematic liquid crystals are fluids of elongated or oblate 
particles, which lack translational order, but exhibit long- 
range orientational orderEla: the molecular axes have a 
common preferred direction, the director. Surfaces and 
interfaces align the nearby molecules and thus favor cer- 
tain director orientations in the bulkcl. This phenomenon, 
known as anchoring, plays Eulsey role in the design of 
liquid-crystal display deviced-M. From a technological 
point of view, one is particularly interested in tailor- 
making surfaces which orient a liquid crystal with an 
arbitrary, well-defined anchoring angle. 

In practice, alignment layers are often produced by 
mechanical rubbing of a polymer coated surface. The 
technique is simple and successful, yet it has drawbacks: 
Low tilt angles between the director and the surface nor- 
mal cannot be achieved easily, and the aljUaament mech- 
anism is still not understood completelyEm. As an al- 
ternative, Halperin and Williams have suggested to use 
swollen brushes of liquid crystalline polymers as align- 
ment layersailj. Their idea was to create a competition 
between the alignment favored by the bare substrate and 
that enforced by the stretching of polymers in a dense 
brush. The interplay between the distortion energy of 
the solvent director field and the conformational entropy 
of the brush was predicted to trigger a second order an- 
choring transition between a phase with planar (homoge- 
neous) alignment at low grafting densities and one with 
tilted alignment at higher grafting densities. The tilt an- 
gle can be tuned by adjusting the grafting density. 

Subsequently, Peng, Johannsmann and Rune have un- 
dertaken first steps -towards an experimental realization 
of such a scenaricO'EJ. A "grafting-from" techniqueE3 al- 
lowed to grow thick side-chain liquid crystalline polymer 
brushes from a substrate covered with surface-attached 
initiators. A competition as required by Halperin and 



Williams can be set up by covering the surface be- 
tween the grafting sites with alkyl chains, which favor 
homeotropic alignment on the bare substrate. The re- 
sulting anchoring behavior has not yet been investigated. 
However, preliminary miscibility studies have shown that 
the liquid crystalline brushes can only be swollen to a 
very limited extent by low molecular weight nematic com- 
pounds, even if the latter are chemically similarEJ. One 
would thus expect that the scenario is altered: The brush 
and the nematic bulk are presumably separated by an 
interface, which will contribute significantly to the an- 
choring properties of the alignment layer. 

This aspect of surface anchoring on grafted liquid crys- 
talline chains is explored in the present paper. We study 
brushes of chains which are too short or too stiff to de- 
velop hairpins and the like, and hence mainly follow the 
director profile. If they attract each other, such short 
chains exhibit collective tihEj. The contribution of the 
chain entropy to the free energy is of minor importance. 
Instead, the elastic energy of director distortions com- 
petes with a tendency of the chain ends at the surface of 
the chain layer to stand up, so that more solvent particles 
can intrude into the interfacial region. As we shall see, 
the interplay between these two factors leads to a novel 
anchoring transition between a phase with tilted align- 
ment and one where the layer of grafted chains aligns the 
bulk vertically (homeotropic alignment), even though the 
chains inside of the layer are tilted. This is demonstrated 
with computer simulations and analyzed within a simple 
theoretical model. 

Our paper is organized as follows: We discuss the sim- 
ulation method and the relevant simulation results in sec- 
tion ||. The theoretical model is introduced and analyzed 
in section HI. We summarize and conclude in section IV. 
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II. MONTE CARLO SIMULATIONS 

We have performed Monte Carlo simulations of a fluid 
of axially symmetric ellipsoidal particles with elongation 
k = cTcnd ond /c 3i do- a idc = 3. The system was confined be- 
tween two hard walls, to which chains of the same parti- 
cles were attached at one end. The grafting points were 
on a regular square lattice. The simulations were con- 
ducted in the NPT ensemble, at constant temperature 
and pressure, and fixed number of particles. Two chain 
monomers and/or solvent particles with orientations u^, 
Uj (|u| = 1), whose centers are separated by a distance 
vector Tij , interact through the pair potential 
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approximates the contact distance in the direction of 
fij|— j= tij I Vjj of two ellipsoids with orientations Uj and 
Ujtj. The anisotropy parameter \ 1S defined as \ — 
(k 2 — 1)/(k 2 + 1). A hard core potential prevents the 
particles from penetrating the walls at z — and z = L z : 
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with d z {6) = cr /2 y/l + cos 2 (0) (n 2 - 1). 

The function d z (9) is the contact distance between the 
surface and an ellipsoid of elongation k and diameter ctq 
oriented with an angle 9 with respect to the surface nor- 
mal. 

Chain monomers are connected by bonds of length 6, 
which are subject to a spring potential with an equilib- 
rium length b — bo and a logarithmic cutoff at | £> — i>o I = 
b,. 



V s (b) = 



o : \b - &o| > b s 



(4) 



The bonds are made stiff and coupled to the orientations 
of the monomers by virtue of a stiffness potential 

V A (0 1 ,6 2 ,9 12 )=k a [4-003(00-008(02)- 2 cos(0 12 )], 

(5) 

which depends on the angles 9\ and 02 between the ori- 
entation of a monomer and the adjacent bonds, and the 



angle 0i2 between the two bonds. The first bond of each 
chain is attached to one of the surfaces. 

We chose the model parameters k s — 10eo/cr 2 , k a = 
10eo, bo = 4.CTo and b s — 0.8tTo. The simulations were 
performed at the temperature T = 0.5eo/^s and pres- 
sure P = 3eo/og. This corresponds to a state well in 
the ncmatic phase: The transition to the isotropic phase 
occurs at the pressure P — 2.ieo/ (J o- ^ ne bulk number 
density was p = 0.313/oq. 

We used a rectangular simulation box of size L|| x Ln x 
L z with periodic boundary conditions in the x and y di- 
rection and fixed boundary conditions in the z direction. 
The Monte Carlo moves included particle displacements, 
rescaling of the simulation box in the z direction (the 
lateral box size L\\ was kept constant in order to main- 
tain a fixed grafting density), and a speciaL-variant of 
configurational biased Monte Carlo movesU : Efl, in which 
chain monomers are turned into solvent particles and new 
chains are grown from the solven t . ., D etails of this algo- 
rithm will be published elsewhereE3llj. 

We have studied systems with roughly 2000 solvent 
particles (the numbers varied slightly in the different 
runs) and up to 242 chains of four monomers in simu- 
lation boxes of lateral size Ln = Via®. The length L z 
of the boxes fluctuated with L z sa 65cto at the highest 
grafting density. The walls orient the solvent particles 
parallel to the surface. We started from an initial config- 
uration where particles all pointed into the x direction. 
This was achieved by applying a strong orienting field 
over roughly 50.000 Monte Carlo steps. After turning 
this field off, the system was equilibrated over at least 1 
million Monte Carlo steps; data were then collected over 
5 million or more Monte Carlo steps. 

The results from these simulations for a range of graft- 
ing densities S will be presented in detail elsewhereEl 
Here, we focus on the anchoring effect that is the subject 
of this paper. Figure shows two configuration snap- 
shots of systems at grafting densities S = 0.34/ctq and 
S = 0.84/oq. The grafted chains are tilted in both cases. 
At the lower grafting density £ = 0.34/ctq, the tilt prop- 
agates into the bulk of the film. At E = 0.84/ctq, the 
orientation of the liquid crystal outside in the bulk is 
perpendicular to the surface, even though the chains still 
retain tilt. 

One might suspect an equilibration problem. In order 
to exclude this possibility, we have prepared an initial 
configuration in which all particles are oriented in the z 
direction. It relaxes into the same structure as Fig. |l|. 

The effect can be characterized more quantitatively by 
inspection of local order parameter profiles. rfThe order 
tensor of a system of n particles is defined bya 



1 ™ q 



(6) 



where Uj denotes the orientation of the particle i as 
above, I the unity matrix, and (g> the dyadic product. 
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FIG. 1. Configuration snapshots (left half of the simu- 
lation box) at grafting densities E = 0.34/<7q and (top) and 
E = 0.84/<j (bottom). Solvent particles are dark, and chain 
particles are bright. 



The largest eigenvalue of this matrix is the nematic order 
parameter S (roughly 0.75 in the bulk of our system), 
and the corresponding eigenvector is the director n. In 
order to obtain profiles of these quantities, we subdivide 
the system in the z direction into slabs of thickness Sz = 
O.lco, and determine the order tensor in each slab. Note 
that the value of S is larger in a slab than in the bulk 
due to the small number of particles in each slab. 

Figs. H a) and b) show averaged profiles of the order 
parameter and of the z-component of the director for the 
same grafting densities as in Fig. |l|. The order parameter 
varies very little throughout the film; it is slightly larger 
in the chain region, due to fact that chain monomers have 
less rotational freedom than solvent particles. The direc- 
tion of alignment varies much more. Close to the surface, 
the particles are aligned parallel to the wall. However, 
the chains reorient the director within a few particle di- 
ameters from the surface. Then follows a region of slower 
changes, where the director gradually stands up. This re- 
gion is restricted to the inside of the chain layer; further 
changes outside of the chain layer are small. In the case of 
S = 0.84/ctq, the last monomers of the chains are almost 
perpendicular to the surface, even though the director is 
tilted inside of the chain region. 

To characterize the chain layer further, Figs. || also 



show the density of chains and the density of chain ends 
(end monomers) . A few chain ends are located very close 
to the substrate at distances z = 1 — 2ao; they belong 
to chains which lie flat on the surface. The other chain 
ends are concentrated at the interface between the chain 
region and the bulk fluid. The higher the grafting density, 
the more efficiently the chain ends are expelled from the 
inside of the chain layer towards the chain-bulk interface. 

Since the cell width L» = 12ao is rather small, quanti- 
tative details of the curves shown in Fig. ^| are presum- 
ably subject to finite size effects. The simulations were 
very time consuming, and we were not able to consider 
different system sizes. In the following, we shall focus 
on the qualitative observation that the chains produce 
homeotropic alignment beyond a certain grafting density, 
even though they themselves remain tilted. 
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FIG. 2. Profiles of the order parameter S (dotted line) 
and the z component (dashed line) of the director n as a 
function of the distance from the surface z for grafting den- 
sities E = 0.34/crg (a) and E = 0.84/crg (b). Also shown 
for comparison is the density profile of chain particles (thin 
dashed line) and the density of chain end particles (thin solid 
line) in units of 1/oq. 
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III. A THEORETICAL MODEL 

In this section, we shall propose a simple theoretical 
model for liquid crystalline chains in a nematic solvent 
which can explain the phenomenon reported in the pre- 
vious section. 

At first sight, the observation that the director ori- 
entation outside of the chain layer differs from that in 
the chain layer is surprising. The "interface" between 
the chain layer and the pure solvent seems to orient the 
particles in a homeotropic way. However, the origin of 
that anchoring force is not obvious. If one admits the 
existence of such an aligning force, one is left with the 
question why the chains are still tilted. In the scenario 
of Halperin and Williamsa the tilt is inherited from the 
bare substrate, which favors planar anchoring. In the 
simulations, the rapid director changes in the vicinity of 
the substrate suggest that the director orientation in the 
chain layer is largely decoupled from that at the sub- 
strate. Alternatively, collective tilt can be induced by 
attractive interactions between chainala. Such interac- 
tions can be mediated by solvent particles due to local 
packing effects, even if the direct interactions are purely 
repulsiveEj. 

Our model is designed to investigate on a simplified 
level the effect of such an interaction. We describe the 
chains by smootb-,differentiable paths R(Z) of length L 
(wormlike chainsEJcJ) , which follow exactly the director 
field of the surrounding solvent: dR/dl = n(R), and are 
attached to the surface at one end. The grafting density 
is measured in terms of a dimensionless area fraction £, 
which is defined as the grafting density times the cross 
section area of one chain. Fluctuations of the director 
are neglected. Moreasier, we do not account for the pos- 
sibility of hairpinstjtj. 

The model has the following four key ingredients: 

1. The elasticity of the nematic fluid. 

2. The translational entropy of the solvent. 

3. The interactions between the chains. 

4. The anchoring energy of the substrate. 

The elastic-energy is given by the Frank free energy 
functionaO~c3 

^ Frank {n(r)} = ^ J dr{K u [V • n] 2 + 

K 22 [n • (V x n)] 2 + K 33 [n x (V x n)] 2 }, (7) 

with the Frank elastic constants K\\ (splay), K 22 (twist) 
and if 33 (bend). 

The translational entropy of the solvent Js calculated 
within the Flory-Huggins approximationOLj 

^ sal ™ t {$(r)} = Ml / dr (1 - $) ln(l - $). (8) 



Here Ub is the Boltzmann factor, T the temperature, vs 
the volume occupied by a solvent particle, and <i>(r) the 
local volume fraction of chain monomers. 

The excluded volume interactions between chains are 
incorporated by the constraint that $(r) must never ex- 
ceed one. The attractive interactions are described by a 
quadratic free energy contribution, 

*W«(r)} = -~ J dr$ 2 , (9) 

where v is an effective interaction strength. Since the 
chains are monodisperse and have no hairpins, their ends 
are located on a smooth, well-defined surface. In the 
vicinity of that surface, the number of interacting neigh- 
bors of a monomer is reduced, which in turn reduces the 
interaction energy. We account for this by adding a cor- 
rection term 

F Burfaoe {$(r);n(r)}=vy ds $ 2 7 ($, n, u s ). (10) 
The integral f ds runs over the surface of the chain 

tJ J surface 

layer, and u s denotes the unit vector in direction of the 
surface normal. The function 7 integrates over the frac- 
tion of "missing" interaction energy in a region close to 
the surface. We take 7 to be proportional to the thickness 
of the affected region, which we calculate as illustrated in 
Figure |[ 7 cx asin(0), where 9 is the angle between the 
monomer orientation n and the surface normal u s , and 
a the distance between two neighbor chains, a oc I/a/^- 
This yields 

7 ($,n,u s ) =7$-V 2 A /T= (n-u) 2 . (11) 

The proportionality constant 7 depends on the form of 
the chain interactions and the actual chain packing at 
the surface, and is presumably of order 1. 




FIG. 3. Surface region of the chain layer where the inter- 
actions between chains are reduced. See text for explanation. 

The free energy contribution ([l0]), ( [Tl"| ) turns out to 
be essential for the effect that we wish to study. On the 
other hand, the considerations which led to eqn. ( J5~l| ) 
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may seem somewhat crude and oversimplified. There- 
fore, we add some comments which justify an ansatz of 
the form ([ll]) on more general grounds: In an anisotropic 
system, one clearly expects the missing neighbor con- 
tribution ([To]) at the surface of the chain layer to be 
anisotropic, i.e., to depend on the angle between the 
surface normal and the director. Since the chains are 
aligned along the director, interchain interactions take 
place in directions perpendicular to the director mainly. 
Hence fewer interactions are missing if the director is 
perpendicular to the surface, and the missing neighbor 
contribution should favor small 0. The ansatz (11), 
•Surface oc sin(0), ensures this in the simplest possible way. 
Other expressions which favor small lead to similar re- 
sults. 

Finally, the orienting effect of the bare substrate is 
described by a quadratic potential 



^ r S ub 3tr ato{n(r)} = W ds(n 

J substrate 



Us 



(12) 



The surface integral J ds runs over the surface of the sub- 
strate. Other substrate-related free energy contributions, 
such as missing neighbor effects of the type ( |l0| ) and sol- 
vent mediated interactions between the chains and the 
substrate are taken to be negligible compared to the an- 
choring energy ([l2|). Furthermore, we neglect the internal 
stiffness of chains, and disregard the possibility that the 
elastic constants may vary locally and depend on the vol- 
ume fraction occupied by the chains. These contributions 
can be incorporated in our model in a relatively straight- 
forward way, but at the expense of having more model 
parameters. The results do not change qualitatively. 

Since fluctuations are neglected, quantities vary only 
in the direction perpendicular to the surface, the z- 
direction, and the director n does not rotate in the xy- 
plane. Inside of the chain layer, the local volume fraction 
<f>(z) of the chains is given by Q(z) — (/n z (z). Outside of 
the chain layer, $(z) is zero and the director is constant. 
Hence the system can be described entirely in terms of 
the profile n z (z) inside of the chain layer. The free energy 
per surface area is given by 



— n z 

^(l-^)ln(l-^)-^) 2 } 

v s n z n z 2 n z J 




The height H of the chain layer is determined implicitly 
by the requirement L = J^ 1 dz/n z (z). 

For practical calculations, it is convenient to perform 
the variable transformation dl = dz/n z (z), i.e., to con- 
sider variations along the chains rather than variations 
in the z-direction. Furthermore, we define the units 
of length and energyta such that (K\\ + Kz%)/2 = 1 



and kBT/v s = 1, and introduce the notation SK — 
(K 33 - Kn)/(K 33 + Kxi). The free energy per sur- 
face area is then expressed as a functional of the profile 
9(1) = arccos(rt z (Z)) of the angle 9 between n and the 
surface normal: 



T{9(1)} 
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and U(x) = (x - 1) ln(l - -) - 

x 2x 

Minimizing this functional yields the Euler-Lagrange 
equation 
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(17) 



and the boundary conditions 



I = 



I = L 



C(9)^- = -2Wcos(9)sm(9) 
dl 

or 9 = arccos(C) (18) 

dl v Vco^(0) 

or = 0. (19) 

Eqn. ( [l7| ) describes the trajectory of a particle in space 
and time I with a position dependent mass C(0), which 
moves in the potential —U(9). The boundary conditions 
pull the particle in the direction of large at I = L and 
in the direction of small at I — 0. Hence the profile of 
decreases monotonically with /. The anchoring angle on 
the chain layer is given by 9(L). 

The properties of the model are most easily discussed 
in the limit of long chains, L — > oo. The two boundaries 
of the chain layer at I = and I = L are then decoupled, 
and there exists a region inside of the chain layer in which 
9(1) is constant, 9(1) = 9. The stationary angle is the 
angle where U (9) has its minimum. It depends on the 
grafting density £ and the chain interaction parameter v. 
If v > 1, the chains exhibits a continuous transition from 
an untilted state at large ( to a tilted state at £ < Cc(v), 
where Cc(v) is determined by the implicit equation 



cm- 



■Kc(«)+ln[l-Cc(«)]}. 



(20) 



In the vicinity of the phase transition, at £ sa Cc(v), 
potential U(9) can be expanded like 



the 
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U{9) = const. - -(C c 



HCc) 9 2 + b(( c ) 9 A (21) 



with 6(C) = (2 - C)/(l - C) + 2/C ln(l - C), and the 
minimum is found at 



- v/2(Cc - C)/Cc 



(22) 



A necessary criterion for the validity of the expansion is 
obviously $ = £/ cos(#) < 1, which is fulfilled if 



C c - C < 2C c arccos 2 (Cc)- 



(23) 



Far from the transition, U{9) takes its minimum close to 
the maximum value of 9 at cos(0)/C 1 and one obtains 



f/(O)-C/(0) = (l-C)[ln(l-C)+</2] 



(24) 



In the case L — > oo, the anchoring angle 6*(L) of the 
nematic film at the surface of the chain layer does not de- 
pend on the anchoring strength W of the substrate. Ex- 
ploiting the integration constant of the Euler-Lagrange 
equation (0), U{8) - ^C(6)(§f = const. = U(5), one 
can rewrite the boundary condition ( fl9| ) at I = L as 

£/(0(L)) - U(9) = g{9) or 9(L) = 0. (25) 
. 2 2 3 (1 + 3/2 tan 2 (#)) 2 



with ff (0) = 2 7 VC 



cos(26») 



U(0)-U(9) 




FIG. 4. Graphical solution of Eqn. (§5|). The an- 
choring angle 6(L) is given by the point where g(9) 
(dashed line) and U(9) — U(6)) (solid line) cross. If 
g(0) is larger than U(0) — U(6), the anchoring angle is 
zero. (The parameters leading to this particular plot were 
C =M.l,v — 1.7,7 — 0.5, SK = 0. Units are defined in the 
textHY 



The problem can be solved graphically as sketched in 
Fig. ||. Since g(9) grows monotonically for all SK < 1, 
the curves g{9) and U{6) — U{9) cross at most at one 
point in the region 9 < 9. If such a point exists, it 
defines the anchoring angle 9(L). Otherwise, the case 
9{L) — applies and the anchoring is homeotropic. The 
transition between homeotropic and tilted anchoring is 
continuous. The grafting density at the transition Q*{v) 
can be calculated by solving U(0) — U(9) — g(0) for a 
given interaction parameter v. Iff is small, C*( w ) is close 
to the tilting transition Cc(^) of the chains, and Eqn. ( p2] ) 
applies. The tilt angle 9 of the chains at the transition 
can be approximated by its asymptotic value at small 
grafting densities C*( v ) 0. Using Eqns. (|2fj| ) and ex- 
panding cos(# ) in powers of 9 , one obtains after some 
algebra 



O) ~ Cc(«)/(1 + 27^/3/(1 + ^)). 



At large v, C* is far from Cc and Eqn. (24) leads to 



C(v) - VT+5K/(2^). 



(26) 



(27) 



Fig. H shows two examples of phase diagrams for the 
parameter values j/y/l + 5K = 0.5 and 0.1. We expect 
that typical values of 7/vl + SK are in that range. At 
high grafting density C and low chain interaction param- 
eter v, the chains stand upright, perpendicular to the 
surface. The thin line marks the continuous transition 
Cc(f) to a phase where they tilt collectively. However, 
the orientation of the director outside of the chain layer 
still remains perpendicular. Only at grafting densities 
below a second critical value C*i v ) < Cc(v) (thick line) 
do the chains induce tilted alignment in the bulk of the 
nematic fluid. 

In agreement with the observations from computer 
simulations reported in the previous section [n[ we thus 
find a region in C — v space where the chains arc tilted, 
but nevertheless align the solvent in a homeotropic way. 
The size of that region increases with 7, i.e., with the 
amplitude of the (anisotropic) missing neighbor effect at 
the surface of the chain layer. It is bounded by two con- 
tinuous phase transitions - one at high grafting density 
to a phase where both the chains and the adjacent fluid 
are untilted, and one at low grafting density to a phase 
where they are both tilted. 

So far we have discussed this phenomenon for the limit 
of long chains, L — + 00. In systems with shorter chains, 
the scenario remains qualitatively similar. If the anchor- 
ing on the bare substrate is planar (W > in Eqn. (fL2|)), 
the sharp tilting transition in the chain region is re- 
placed by a smoother crossover from distiiictly tilted to 
roughly perpendicular chain orientationO, and the an- 
choring transition shifts to higher grafting densities. One 
still finds a region where markedly tilted chains align the 
nematic solvent in a homeotropic way. 
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0.1 (b) 



( space (chain interac- 
0.5 (a) and 



7/Vl + SK = 0.1 (b). Thin solid lines indicates tilting tran- 
sition in the chain region, and thick solid line the anchoring 
transition in the nematic fluid. Short dashed line indicates 
the approximation ( ^ ) and long dashed line the approxima- 
tion ( p7| ) (outside p£ the frame in the case of (b)). Units are 
defined in the textLj. 



IV. SUMMARY AND DISCUSSION 

To summarize, we have described a new type of an- 
choring transition on grafted layers of main-chain liquid- 
crystalline chain molecules. It is driven by the compe- 
tition between attractive chain interactions, the transla- 
tional entropy of the nematic solvent in a region close to 
the surface of the chain layer, and the elastic energy of 
the nematic fluid which opposes rapid director changes. 
Inside of the chain layer, the chain interactions make the 
chains tilt collectively. Close to the surface, the inter- 
actions are reduced, and an increased number of solvent 
particles enters and swells the chains. As a result, chain 
ends in the surface region stand up. The anchoring an- 
gle imposed by the chains on the adjoining nematic fluid 



is thus smaller than the actual tilt angle of the chains. 
Above a critical grafting density £*, it drops to zero, de- 
spite the fact that the chains are still tilted. We have 
observed this effect in computer simulations of a model 
system with swollen grafted chains, and rationalized it in 
terms of a simple theoretical model. 

The new anchoring transition is controlled by the graft- 
ing density of the chains. In that sense, it i*. similar 
to that discussed by Halperin and Williamsljil^l. How- 
ever, the phase transition and the underlying mechanism 
is very different. In the case studied by Halperin and 
Williams, the chains are long and flexible enough to sup- 
port many hairpins. The solvent is taken to be very good, 
i. e., chain monomers effectively repel each other. The 
transition is driven by an interplay between the confor- 
mational entropy of the chains, the anchoring force of 
the substrate on the solvent, and the elastic energy in 
the brush. The transition connects a phase with planar 
anchoring and one with tilted anchoring. 

In the present paper, we consider a different regime: 
The chains are shorter or stiffer, they have no hairpins, 
and their conformational entropy is negligible. The sol- 
vent is not good, chains effectively attract each other and 
there exists a well-defined interface between the chain 
layer and the solvent. The anchoring behavior is con- 
trolled by the structure of that interface. The influence 
of the bare substrate is minor. The transition connects a 
phase with tilted anchoring and one with perpendicular 
anchoring. 

One might ask whether there could be a second tran- 
sition between planar and tilted anchoring, at a lower 
grafting density, which would resemble that discussed by 
Halperin and Williams. Our theoretical model does not 
predict such a transition. In the simulations, the situa- 
tion is more complicated. The data indicate the presence 
of a first order transition, which has nothing to do with 
hairpins, but is still related to the organization of the 
chains. ThiB. will be described in detail in a separate 
publicatiorMN. 

Experimentally, one finds that liquid crystalline poly.- 
mers are not swollen very well by nematic solventsEJ. 
We thus believe that our effect can be observed in real 
systems. The grafting density must be high to achieve 
low anchoring angles 9, on the other hand, the chains 
need not be much longer than their persistence length. 
Moreover, our results indicate that the grafting density 
required for the transition decreases with decreasing sol- 
vent quality. Brushes with high grafting densities can 
be prepared with the "grafting from" technique devel- 
oped by Riihe et aLJ. This technique has been used in 
the experiments of Peng et altlrll3 mentioned in the in- 
troduction. However, Peng et al studied brushes with 
liquid crystalline side chains, where the mesogenic units 
are preferably oriented perpendicular to the backbone of 
the chains. The swelling of the brush in the surface re- 
gion thus promotes planar anchoring, and no transition 
to homeotropic anchoring can be expected. The effect 
discussed in the present paper will presumably transpire 
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in a similar experiment with short main-chain liquid- 
crystalline polymers. 

One necessary requirement is of course that the ori- 
entation of the chains is homogeapous. This turned out 
to be a problem experimentallyliil, the brushes tend to 
exhibit planar multidomains. However, recent work indi- 
cates that an appropriate treatment of the bare substrate 
before growing the brush can force the chains into one 
monodomainllj. If brushes could be grown with graft- 
ing densities close to the transition density ^* , it would 
be possible to prepare surfaces which anchor the solvent 
at any given small tilt angle 8*. This would potentially 
have applications in the design of liquid crystal display 
devices. 

Future theoretical work will have to study in more de- 
tail the chain length dependence. The theory presented 
here has been devised for long chains - for example, 
they have been described by homogeneous strings - but 
the chains in the simulations were very short with only 
four monomer units. Chains in real brushes are usu- 
ally polydisperse. The influence of polydispersity on the 
anchoring behavior is not clear. Thus far the analyti- 
cal model calculations have disregarded the possibility 
of lateral structure and lateral fluctuations. These will 
certainly become important for polydisperse brushes and 
more generally at low grafting densities. 
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